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Abstract 



It is shown analytically that the Dirac equation has no normalizable, time-periodic 
solutions in a Reissner-Nordstrom black hole background; in particular, there are no 
static solutions of the Dirac equation in such a background metric. The physical 
interpretation is that Dirac particles can either disappear into the black hole or escape 
to infinity, but they cannot stay on a periodic orbit around the black hole. 

1 Introduction 

In recent years, there has been much interest in the gravitational collapse of matter to 
a black hole. Although both analytical and intensive numerical studies (see e.g. ]|]) 
have given some understanding of how the event horizon and the singularity form, little is 
known about the asymptotic form of the black hole as t — > oo. This is mainly due to the 
fact that standard numerical methods become unreliable after the solutions have formed 
singularities. Since all matter on a microscopic level is formed out of Dirac particles, 
it seems especially interesting to study the asymptotic collapse of a "cloud" of spin-|- 
particles. As a first step towards this goal, in this paper we study Dirac particles in a 
Reissner-Nordstrom background field. 

We remark that considerable work has been done in the study of quantum mechanical 
wave equations in the presence of black holes. The papers which are most related are 
[||, Q, where a massless Dirac particle is considered in a Schwarzschild metric background, 
and asymptotic completeness is shown for the scattering states near the event horizon and 
at infinity. However, the most physically interesting case of a massive Dirac particle near 
a charged black hole has not yet been considered. As we will see here, both the rest mass 
of the Dirac particle and the charge of the black hole lead to interesting physical effects 
and require new analytical tools. 

In polar coordinates (t,r,$,ip), the Reissner-Nordstrom metric has the form 
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where q is the charge of the black hole and p its (ADM) mass. Furthermore, we have an 
external electromagnetic potential A of the form A = (—</>, 0) with the Coulomb potential 



4>(r) 



(1.2) 



r 



If q < p, the metric has two horizons; this is the so-called non-extreme case. If q = p, 
the metric has only one horizon at r = p; this extreme case describes a black hole at zero 
temperature; cf. || ||, [?]]. For q > p, the metric does not describe a black hole, and thus 
this case will not be considered. 

We describe the Dirac particles with "classical" wave functions (i.e. without second 
quantization). Both the gravitational and electric fields are coupled to the Dirac particles. 
We do not assume any spatial symmetry on the wave functions. Near a collapsing black 
hole, it seems reasonable that that some of the Dirac particles could get into static or 
time-periodic states. Our main result is to show that this is not possible. 

In the following we will restrict to time-periodic solutions, noting that static solutions 
are a special case. For classical point particles, the time-periodic solutions describe closed 
orbits of particles rotating around the black hole. Our goal is to investigate how this 
classical picture changes by the introduction of relativistic wave mechanics and spin. Since 
the phase of the Dirac wave function \& is of no physical significance, we say that ^ is 
periodic with period T if 



for some real J7. Our main result in the non-extreme case is the following theorem: 

Theorem 1.1 In a non-extreme Reissner-Nordstrom black-hole background, there are no 
normalizable, periodic solutions of the Dirac equation. 

In the extreme case, we prove a slightly weaker statement: 

Theorem 1.2 In an extreme Reissner-Nordstrom background, every normalizable, time- 
periodic solution of the Dirac equation vanishes identically for r > p. 

This surprising result shows that the classical picture breaks down completely; for Dirac 
particles, there are no periodic solutions. This means that Dirac particles which are 
attracted by a Reissner-Nordstrom black hole either "fall into" the singularity or escape 
to infinity, but they cannot stay on a periodic orbit around the black hole. The result can 
also be applied to the Dirac particles of the matter in the gravitational collapse; it then 
indicates that all the matter must eventually disappear in the black hole. 

Basically, our result is a consequence of the Heisenberg Uncertainty Principle and of 
the particular form of the Dirac current. As a preparatory step, we analyze the behavior of 
the Dirac wave functions near the event horizon and we derive conditions which relate the 
wave function outside and inside the horizon. It is essential for our methods and results 
that the particles have spin. This shows that the spin is an important effect to be taken 
into account in the study of gravitational collapse. 

In the remainder of this section, we give some basic formulas needed to describe Dirac 
particles in curved space-time (for a more detailed introduction to the classical Dirac 
theory in curved space-time see ||]). In this paper, the Dirac equation is always of the 



(1.3) 



form 




(1.4) 
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where m is the rest mass of the particle, A = Ajdx J the electromagnetic potential, and e 
the electromagnetic coupling constant (see |J, [10] for a derivation of this equation). The 



Dirac matrices G 3 (x) are real linear combinations of the usual 7-matrices. We work in the 
Dirac representation 

M° -0 ■ M-°*< 0) ■ <=™ <-) 

where a % denote the Pauli matrices. The Dirac matrices are related to the Lorentzian 
metric via the anti-commutation relations 

9 3k {x) = \{<?(x), G k (x)} • (1.6) 

The term VjG J in ( |1.4p is the divergence with respect to the Levi-Civita connection; it 
can be easily computed via the standard formula 

V S (P = -j=dj(J\dG>) . (1.7) 



For the normalization of the wave functions, one takes a space-like hypersurface TL with 
normal vector field v and considers the scalar product 



(tf|$)= / i&G 3 $v j dn , (1.8) 
JH 

where ^ = ^*7° is the adjoint spinor, and where dp, is the invariant measure on TL induced 
by the Lorentzian metric. On solutions of the Dirac equation, we impose the normalization 
condition 

(* I *) = 1 . 

Current conservation 

V 3 fG j f = (1.9) 

implies that this normalization condition remains unchanged if the hypersurface TL is 
continuously deformed. 

2 The Dirac Operator in a Schwarzschild Background 

We begin by analyzing the Dirac operator in a Schwarzschild background metric. Our 
aim is to analyze the behavior of the spinors near the event horizon. To do this, we must 
consider the Dirac equation in different coordinate systems. 

2.1 The Dirac Operator in Polar Coordinates 

In polar coordinates (t, r, $, ip), the Schwarzschild metric is 



ds 2 = \1 - y J dt 2 - M - j dr 2 - r 2 {d$ 2 + sin 2 ■& dip 

where p is the (ADM) mass. The metric has an event horizon at r = 2p. In order to derive 
the Dirac operator, we first choose Dirac matrices G 3 {x) satisfying the anti-commutation 
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relations ( |1.6[) . The Dirac operator is then obtained by calculating the divergence (|1.7| ) 
and substituting into ( p.-4| )P[ 

Outside the horizon, we can satisfy the anti-commutation relations (|1.6| ) by choosing 
the Dirac matrices in the form 

G* = i 7 * , G r = Sf , G^ = ^ , G^ = ^ {r>2p) (2.1) 

with 

S(r) = l-^ 1 
r 

where 7*, j r , 7 , and 7^ are the "7-matrices in polar coordinates" 

,0 

cost? + 7 1 sint? cosy? + sin$ s'mip (2.2) 

(2.3) 

(— 7 1 sin 93 + 7 2 cos 99) . (2-4) 

The divergence of the Dirac matrices is computed to be 

Vj<? = (s' + 2 - (5 - 1)) 7 r • 
Substituting into (|1.4j) , we obtain for the Dirac operator, G ou t, in the region r > 2p 

e - " I y I + ^ ( ,s I + 7 (s - " + I s ') + ^ + "< e h ■ < 2 - 5 > 

For the normalization, we integrate over the hypersurface t = const; i.e. 

(*|*)^:= I (Wm,*) S- 1 d 3 x , (2.6) 
JtR 3 \B 2p 

where F>2 P denotes the ball of radius 2p around the origin. This normalization integral is 
problematic near the event horizon, as will be discussed in detail later. Inside the horizon, 
we must take into account that the radial direction r is time-like, whereas t is a space 
coordinate. So in this region, to obtain the Dirac matrices, we reverse the roles of the 
matrices 7* and 7**, 

G t = ^Y , G r = -S-/ 1 , G & = -f ,t} , G* = ^ (r<2p). (2.7) 
The divergence of the Dirac matrices now has the form 



VjG j = --f - (s' + -s)^ 



r 



We point out that the choice of the Dira c m atrices is not canonical; there are different real linear 
combinations of the 7-matrices which satisfy (|l.6|) . But the Dirac operators corresponding to different 
choices of the Dirac matrices are equivalent in the sense that they can be obtained from each other by a 
suitable local transformation of the spinors (see e.g. Q). For this reason, we can simply choose the G J in 
the way which is most convenient to us. 
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Thus the Dirac operator, Gi n , in the region r < 2/3 is given by 



According to ( |1,8| ), the naive extension of the normalization integral ( |2.6| ) to the interior 
of the horizon is 

(* | vp)f n := / (*Ym,£) S' 1 d s x ; (2.9) 

this will also be discussed in detail later. 

Notice that as a particular freedom in the choice of the Dirac matrices, the signs in 



( |2.lD and (2^) are arbitrary. As remarked above, this arbitrariness can be compensated by 
a suitable local transformation of the spinors. However, this transformation of the spinors 
may change the sign of the scalar product ( |1.8|) . This is a subtle point which needs some 
explanation. Assume that we consider the space-like hypersurface outside the horizon 

Hi = {t = const, r > 2p} . (2.10) 

Its normal vector field v is only determined up to a sign. Depending on whether we choose 



for v the future or past directed normals, the corresponding scalar product (1.8) will (for 
a fixed choice of the Dirac matrices G 3 ) be either positive or negative (semi-)definite. 
However, the overall sign of the scalar product is of no physical relevance; e.g., we could 



just redefine (1.8) by inserting a minus sign. In order to fix the sign convention, we will in 
the following always assume that the scalar product (|1.8j ) is positive for the future-directed 
normal vector field (this convention is consistent with our choices ( |2.1| ) and ( |2.6[) ). The 
situation becomes more interesting if we also look at the region inside the horizon. For 
this, we consider the "cylindric" space-like hypersurface 

n 2 = {r = r , i < i < *i} (2.11) 
for some fixed ro < 2p and to < t\. A short computation shows that, for our choice of 



the Dirac matrices (|2.7| ), the scalar product (|l.8|) corresponding to H2 is positive if we 
choose for v the inner normal (pointing towards the singularity at r = 0). According 
to our sign convention, this means that the inward radial direction points to the future. 
Thus the particles "fall into" the singularity as time progresses, and we have a black hole. 
On the other hand, we could have chosen the Dirac matrices such that the scalar product 
corresponding to Tii is positive for the outer normal (e.g. by changing the sign of G r in 
(p-Tp). In this case, increasing r would correspond to going forward in time, and we would 



have a white hole. Notice that this argument is consistent with time reversals. Namely, 



the replacement t — > — t forces us to change the sign of the scalar product (1.8) (in order 



that (|l.8|) is still positive for H = Hi and future-directed normals). As a consequence, the 
scalar product corresponding to H.2 changes sign. This means that black holes become 
white holes and vice versa. We conclude that the Dirac operators G ou t and Gj n distinguish 
between a black and a white hole. This is a peculiar effect of the Dirac operator. It is quite 
different from e.g. the wave operator describing scalar fields (the Klein-Gordon operator), 
which does not determine the direction of time inside the horizon. 

Our description of the spinors in polar coordinates is not quite satisfactory. First of all, 
the normalization integral inside the horizon, (|2.9| ), is not definite. This is a consequence 
of the fact that the t- variable is space-like inside the horizon. From the mathematical 
point of view, this is no problem; it seems tempting to just integrate across the horizon by 
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Figure 1: Kruskal Coordinates 



adding ( |2.6| ) and ( p.9[ ). On the other hand, it is a conceptual difficulty that the integrand in 
( |2.9D is not positive and therefore does not have the interpretation as a probability density. 
Furthermore, the Dirac equations corresponding to G ou t and G in separately describe the 
wave functions outside and inside the horizon. But it is not clear how to match the wave 
functions on the horizon. For a better understanding of these issues, it is useful to remove 
the singularity of the metric on the horizon by transforming to Kruskal coordinates. 



2.2 Kruskal Coordinates 

According to [pLlJl , we introduce Kruskal coordinates u and v by 



1 e 4 P C osh ( — ) for r > 2p 

u = { . Yfl (2-12) 



f r ft 

1 e^p sinh I — I for r <2p 

2p \Ap ' 

~T r ft 

1 e^p sinh I — I for r > 2p 

2p V4p ' 

T r ft 

1 eM cosh I — ) for r < 2p 



(2.13) 



2p \4 P/ 

The regions r > 2p outside and r < 2p inside the horizon are mapped into 

0\ = {u> 0, |v| < u} 

and 

h = {v > 0, \u\ < v, v 2 - u 2 < 1} , 

respectively (see Figure |l|). The horizon r = 2p corresponds to the origin u = = v, 
and the lines v = ±u are reached in the limit t — > ±oo. Finally, the singularity at r = 
corresponds to the hyperbola v 2 — u 2 = 1, v > 0. 
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In Kruskal coordinates (v,u, 1?, <p), the Schwarzschild metric takes the form 
ds 2 = f~ 2 {dv 2 - du 2 ) - r 2 {d$ 2 + sin 2 dip 2 ) with 

r 

This metric is regular except at the singularity v 2 — u 2 = 1; it can be extended to the 
entire region v 2 — u 2 < 1. 

Since the metric is regular at the origin, we can smoothly extend the Dirac operator 
across the horizon. To do this, we simply view v and u as the time and space variables, 
respectively. We choose for the Dirac matrices 

G v = / 7 * , G u = f 7 r , = ^ , G v = 7^ . 

A straightforward computation yields for the Dirac operator 

g = y(/i|- + i/( ft ,r)-ifl B /)+7 p (/i^ + J(/(aur)-i)-ia I ,/) 

+ i7% + »7 1p 3 v • (2-14) 

The normalization integrals ( |2.6D and ( |2.9D on the surface t = correspond in Kruskal 
coordinates to the integral ( p..8|) with 

H = { u = 0, < v < 1} U {v = 0, u > 0} . 

We choose the normal v as in Figure Using the current conservation (|1.9| ), one can 
continuously deform the hypersurface TL without changing the value of the normalization 
integral. In particular, we can avoid integrating across the horizon by choosing the hy- 
persurface Ti in Figure |. This is a major advantage of Kruskal coordinates; it gives a 
physically reasonable positive normalization integral even inside the event horizon. How- 
ever, this method must be done with care when the considered solution of the Dirac 
equation has singularities near the origin. Unfortunately, our time-periodic solutions of 
the Dirac equation will, after transforming to Kruskal coordinates, in general be highly 
singular at the origin. Therefore the deformation of the hypersurface as in Figure |l| would 
be problematic, and we will not use this method. In order to avoid any difficulties of the 
normalization integral near the horizon, we shall only consider the normalization integral 
outside and away from the event horizon. 



2.3 Transformation of the Dirac Operator 



We now consider how the Dirac operator ( |2.5[ ), Q2.8| ) in polar coordinates transforms 
into the Dirac operator ( |2,14j ) in Kruskal coordinates. This transformation consists of 
transforming both the space-time coordinates and the spinors. For clarity, we perform 
these transformations in two separate steps. Under the transformation of the space-time 
coordinates, the partial derivatives transform as 

1/5 5 
— u tt - + v t; — 
4p V dv du, 

1/5 5 , 

v — + u — I for r > 2p 





dt 


dv d 
dt dv 


+ 


du 
~di 


5 
du 


d 


dv d 


+ 


du 


d 


Or 


dr dv 


dr 


du 



4p S 2 V dv du, 

1 ( 9 5 . 

v — — + u — — tor r < 2p 



4p S 2 V dv du 
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Substituting into (2.5) and (2.8) gives for the Dirac operators G out 
coordinates 



and Gi n i n Kruskal 



G 



Gi, 



(u7* + wy r ) — 



du 







(U7 + U7 ) — 



-S +-S' 
r 2 



- 7 + *7 + *7 V ^~ 



(2.15) 



(2.16) 



These Dirac operators do not coincide with ( |2.14 ), and we must therefore perform a 
further transformation; namely a transformation of the spinors. Under general coordinate 
transformations, the wave functions transform according to 



^(x) -> U{x) ^>{x) ; 
in the case considered here, U is the time-dependent (4 x 4) matrix 



U(t) = cosh 



1 + sinh 



7*7 r 



Under this transformation, the Dirac operators G, 



out/in 



t 

transform as 



G, 



out /in 



^ ^out/in^ 



-1 



(2.17) 



(2.18) 



(2.19) 



This gives the Dirac operator ( |2.14| ) in Kruskal coordinates, 

G = U I/" 1 = 17 G in U- 1 



and this can be verified as follows: Under the transformation ( 2.19 ), the Dirac matrices 
behave like 

G j (x) -» U(x)G j (x)U(x)- 1 . 

Using the form of the Dirac matrices in ( |2.15[ ), ( p,16| ) and the explicit formula ( 2.18; ), a 
short calculation shows that the Dirac matrices of the operators £7G out /j n i7 _1 coincide 
with the Dirac matrices /7*, fj r , 7 , and 7^ in ( |2,14| ). According to (1.4), the Dirac 
operator in the gravitational field is formed from the Dirac matrices and their covariant 
derivatives; it is thus completely determined by the Dirac matrices. Therefore, the op- 
erators t7G out / in f/ _1 must coincide with G. (One can also verify explicitly that the zero 
order terms of the operators f7G out / in f/ _1 and G are equal. This is a longer computation, 
however.) We conclude that the Dirac operators G ou t and Gj n can be identified with the 
Dirac operator G in the region 0\ U I\ = {u + v > 0, v 2 — u 2 < 1}. 

We remark that it is not possible to map the interior of the horizon into the region 

h = {v < 0, |«| < -v, v 2 - u 2 < 1} 

and still match the Dirac operator G; n with G, because this would contradict the fact 
discussed in the previous section that the Dirac operator distinguishes between black and 
white holes. 

Finally, we note that the transformation ( 2.17[ ),( p.l9| ) of the spinors can be viewed as 
a local C/(2, 2) gauge transformation; see ||. 
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2.4 Matching of the Spinors on the Horizon 

We now come to the question of how the wave functions inside and outside the horizon are 
related to each other. For this, we analyze the behavior of solutions of the Dirac equation 
at the origin in Kruskal coordinates. After transforming back to polar coordinates, this 
will give matching conditions for the wave functions on the event horizon. The physical 
situation which we have in mind is a Dirac particle attracted by a Schwarzschild black 
hole. It suffices to do the matching for static solutions, (and not time-periodic solutions), 
since in Section we reduce the problem to static solutions. 
Let f be a static wave function, i.e. in polar coordinates 

We assume that \I/(r, $,(p) is a smooth function both inside and outside the horizon; 
i.e. in the regions r < 2p and r > 2p; this will be justified later by a separation of 
variables technique. Furthermore, we assume that ^ is a solution of the Dirac equations 
(Gi n — m)^ = and (G out — m)^ = 0, respectively. According to the transformation rules 
fl2.17l ), (|2.18|) , the wave function \& in Kruskal coordinates takes the form 

V(u,v,$,tp) = U(t) e~ iuit V(r,0,(p) , (2.20) 

where r and t are given implicitly in terms of u and v by 

— -l] e a? = u 2 -v 2 (2.21) 
2p J 



t 



4p 



tanhf — ) = { I (2.22) 



for r > 2p 



u V 

— tor r < 2p . 

v 



With this procedure, is only defined in 0\ U lj, the upper right half of the Kruskal 
domain; it solves the Dirac equation 

(G-m)t> = (2.23) 

in the open set 0\ U I\. If ^ is to be a physically reasonable solution of the Dirac equation, 
it must be possible to extend it to the entire region v 2 — u 2 < 1 between the two hyperbolas. 
If this extended wave function was not zero in the region u + v < 0, our system would be 
connected to a white hole or to another universe (through a worm hole), and the Dirac 
particle would have a certain probability to be in these extensions of space-time. Since we 
are only interested in black holes, this is not the situation we want to consider. Therefore 
we demand that the extension of must vanish identically in the half plane u + v < 0. We 
conclude that in Kruskal coordinates, we must analyze a solution ^ of the Dirac equation 
fl2~23l) of the form 

f U{t) e"** *(r, *,<p) for u + v > 0, u + v 
v ' ' ' ^' | for u + v < . v ' 

This wave function may be singular on the lines u = ±v; in this case, ^ must solve the 
Dirac equation in a generalized weak sense. 

For the calculation of the weak derivatives of ^, we rewrite the wave function in the 
form 

* = 6(u + v) 0(u - v) $>o + ®(v + u) Q(v - u) , 
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where = ^\Oi an d ^1 = are the components of ^ outside resp. inside the horizon 
(0 denotes the Heaviside function O(x) = 1 for x > and 0(x) = otherwise). Since VP 
satisfies the Dirac equation in 0\ U Ji, we need only consider the singular contributions 
on the lines u = ±v. A formal calculation gives 

= (G - m) * 

(2 = 4) fi ( 7 * + Y) S(u + v) (6(u - u) * + 0(u (2.25) 
-/i ( 7 * - 7 r ) <5(u - u) 0(« + v) (l-o - 17) • (2.26) 

If and VP/ were smooth up to the boundary of 0\ resp. Ii, this equation would be 
well-defined in the distributional sense. In general, however, $>o and VP/ might be sin- 
gular in the limit u — > ±v. In order to treat this general case, we multiply (Bl . fB) 



with test functions 77(71, v) which, as u — > ±v, decay so fast that the integral over the 
resulting expression is well-defined. Since the matrices (7* + 7 r ) and (7* — 7 r ) are linearly 
independent, we get the two conditions 

77 / S(u + v) (7* + 7 r ) (0(u -v)V + ©(« - u) */) du cfo = (2.27) 

IR 2 

/ r ? /5(n-w)(7*-7 r )0(u + u)(^ o -1/)^(iv = . (2.28) 

JIR 2 

In fl2.27p , we compensate the possible divergence of VP for u — > —v by choosing 77 in the 
region Oi U I\ to be of the form 77101U/1 = (1 + 1(7* + l r )^\) 9 with a smooth function 



g. Then the integrand in (2.27) is of the form 5(u + v ) x g x (bounded function), and 
the integral makes sense. Since g is arbitrary, we conclude that the integrand on the line 
u = —v must vanish, which implies that 

lim (7* + 7 r ) V(u,v,d,<p) = . (2.29) 

U — > — V 

In fl2.28| ), we can proceed similarly; namely, if VP is singular on the line u = v, we com- 
pensate the divergence of the integrand in ( |2.28 ) by choosing 77 to have an appropriately 



fast decay near the line u = v. One must keep in mind, however, that 77 cannot be chosen 
independently in 0\ and I±, because the smoothness of 77 on the line u = v may impose 
restrictions on 77. For example, if \Po an d VP/ have poles near u = v, 

u + £, i?, tp) = ci(u,ti,<p) e" p + • • • 
^o(u,u - e,$,ip) = c 2 {v,,$,tp) e~ q + ■ ■ • (e > 0), 

then we must choose 77 in the form 

7?(n, u + e, 0, tp) = c 3 (u, t?, tp) e max (P-5) + . . . (e > or e < 0). 

Thus the asymptotic behavior of 77 near u = v in 0\ and I\ must be the same. In the 



integral ( [2.28D , this means that the leading order singularities of VPo and VP/ may cancel 
each other for any choice of 77. Therefore, the condition for the leading order singularity 
takes the form 

(7* - 7 r ) (VP(u, tj + e, 7?, <^) - VP(u, u - e, 1?, 99)) 
= o(l + |(7* -7 r ) ^(7j,7j + e,7?,yj)|) ase^O. (2.30) 
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If the singularity of \& on the line u = v is worse than polynomial or of different form, 
there may be no obstructions for the choice of rj in 0\ and 1%. In this case, ( |2.30| ) will still 
be a necessary condition. It will no longer be the strongest possible condition, but this is 
irrelevant for our purposes. For simplicity, we will use ( |2.30| ) in the general case. 

Next we evaluate the conditions ( |2.2S| ) and ( |2.30D for our wave function ^(u, v, (p) 
in ( gjg) . Using ( pl^) , we have in 0\ U h 



( 7 * - 7 r ) V(u,v,0,<p) 



(7* + 7 r )e^ e~ iujt $>(r,$,(p) 
(7* 



7 r ) e - *? e~ iujt V(r,d,<p) 



The explicit formulas fl2.12|) and fl2.13| ) enable us to write the time-exponential in terms 
of u and v as 

u + v 



±-£- 

e g p 



±1 



Using the relation (|2.2l| ) between r, u, and u, the condition ( 2.29f) in polar coordinates 
takes the form 

(2.31) 



lim(7* +7 r ) \e\*^(t,2p + £,■&, <p) 







Similarly, condition ( 2.30| ) can be written in polar coordinates as 

( 7 *-7 r )N^ (*(2p + £,#,</>) - ¥(2p-e,<? >¥ >)) 
= o(l + |e|3 |(7*-7 r )^(2/) + e,T?,v9)|) as e -> 0. 

In order to simplify this formula, we consider the decomposition of \& in the form 
| e |3 $ = 7 * ^y + y) | e |i ^ + ( 7 * - 7 r ) | £ |i ^ 



(2.32) 



(2.33) 

According to condition fl2.31| ), the first summand in the bracket in (2.33) vanishes on the 
horizon r = 2p. Since the matrix 7* is invertible, we conclude that ^ and (7* — 7 r ) v I' are 
of the same order on the horizon. Thus we can take out the matrices (7* — 7 r ) in ( |2.32j ) 
and finally obtain the equivalent condition 

|e|* (tf(t,2p + e,i? >¥ >) - tf(t, 2p - e, 0, p)) = o(l + |e|i ¥(t, 2p + e, 0, 99)) (2.34) 

as e — > 0. The relations (|2.3l| ) and ( |2,34j ) are our matching conditions. 

We briefly explain what these matching conditions mean, without being mathemati- 
cally precise. First of all, we point out that the matrix (7* + 7 r ) in the first matching 
condition (2.31) is not invertible. Therefore, ( [2.31 ) does not imply that |e|* ^f(2p + s, $, ip) 
goes to zero in the limit e — > 0; in general, this limit need not even exist. Although the 
matching conditions have a quite special form, they can be understood intuitively if one 
considers the Dirac current in polar coordinates. We first look at the total normalization 
integral M + CT: 



(* I *)out + I *)in 

*(7 r + 7*)* dp 

B 2p 



/ ^7** dp + / *7*^ dp . 

Jb 2p Jm 3 \B 2p 



The condition ( pT3l| ) ensures that the integral of the first summand is small near the 
horizon. Using the matching condition ( |2.34 ), one sees that the integrals in the second 
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and last summands behave similarly near the horizon. Because of the opposite sign of the 
second and third summands, this tends to make the normalization integral finite even if 
\I/ is singular on the horizon (if the current had a pole, for example, one could define the 
normalization integral as a principal value). Thus our matching conditions "regularize" 
the normalization integral across the horizon. It is also interesting to look at the current in 
radial direction. For this, we consider the normalization integral through the hypersurface 
7^2; (2.11). For the outer normal u, this gives inside the horizon 



(y\V) n . 2 = - [ tyrf^dn (r<2p). (2.35) 



I Hi 

For r > 2p, on the other hand, we get the expression 



(* | ^)h 2 = \ *7 r * dp, 

= ( *(7 r + 7*)*d/x - / *7**^ {r>2p). (2.36) 
J Hi JHi 



According to ( 2.31| ), the first integral in fl2.36j ) is small near the horizon r = 2p. The 



matching condition Q2.34 ) gives that the second summand in ( |2.36| ) behaves similar to 



( |2,35| ) near the horizon. Thus our matching conditions tend to make the normalization 
integral through Ti.2 a continuous function in ro on the horizon r$ = 2p. Since the in- 
tegrand of the normalization integral has the interpretation as the "probability density" 
or "probability current," this means physically that a particle which disappears in the 
event horizon must reappear in the interior of the horizon. This is in accordance with our 
physical assumption that there are no other universes or white holes where the particle 
could disappear into or emerge from. 



3 Separation of the Angular and Time Dependence 

We next study Dirac particles in the external Reissner-Nordstrom background fields ( |l.lf ), 
( |1.2[ ). Since the external fields are spherically symmetric and time-independent, we can 
separate out the angular and time dependence of the wave functions via spherical harmon- 
ics and plane waves, respectively. This is done in a manner similar to the central force 
problem in Minkowski space (see e.g. ||l2|| ). 

We start with a compilation of some formulas involving the angular momentum oper- 
ator L = —i(x x V) (see e.g. [l3|]). Its square is 

L 2 = -A S 2 = L + L- + L\ - L z = L^L + + L 2 Z + L Z 

with L± = L x ± iL y . The spherical harmonics Yj fc , I = 0, 1, . . ., k = —I, . . . ,1 are simulta- 
neous eigenfunctions of L 2 and L z , namely 

L 2 Yf = l{l + l)Y t k , L z Y t k = kY t k . (3.1) 

They are orthonormal, 

/ Yt*Y$ = 5 lv 5 kk ' , 
Js 2 

and form a basis of L 2 (S 2 ). The operators L± serve as "ladder operators", in the sense 
that 

L±Y l k = + 1) - k{k ± 1) Y k±1 . (3.2) 
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In preparation for the four-component Dirac spinors, we consider two-component Pauli 
spinors. In analogy to ( |2.2[ )-(EO|), we denote the "Pauli matrices in polar coordinates" by 
cr r , a , and a v ; i.e., 

a r = a 3 cos i? + a 1 sin $ cos <p + a 2 sin $ sin ip 

a® = - (—cr 3 sini? + a 1 cost? cos ip + a 2 cost? sin ip 

1 

r sini? 



a v = : — - f— a 1 sin ip + cr 2 cosy? 

- sin 1/ V 



We have 



and thus 



0*84 + 0*8^ = aV - a r d r = — (ax)(aV - a r d r ) 

T V 

= —{rd r + ia(xx V) - rfl^) = -— aL , (3.3) 

7" V ' T 

at = -r a r (<r% + t^SJ . (3.4) 



For i = 5> |> ■ ■ ■ and & = —J, — j + 1, . . . , j, we introduce the two-spinors 



j + k k-i ( 1 \ j - k fc+± ( 



2 



2j j-3 v / y 2j j-3 v 1 



fe / j + i-fe fc-§ / 1 \ _ j + i + k k +i / \ 
^■+1 " V 2i + 2 i+ 5 V / V 2 i + 2 V 1 / ' 

These spinors form an orthonormal basis of L 2 (S 2 ) 2 . They are eigenvectors of the operator 
K = aL + 1. More precisely, (|3.1| ) and ( |3.2| ) imply that 

K *U - ("C *U = « + l)xU (3.5) 

Kx k j+ l = -U + \)x] +i ■ (3.6) 
Furthermore, multiplication with a r again gives an eigenvector of K; namely 



^ r - r <r p (a* a r % + cx^ a r d v ) 

-'2 J a 

J 2 J 2 



Taking into account the normalization factors, we obtain the simple formula 

° r *U = *U ■ (3 - 7) 

Finally, we choose for the Dirac wave functions the two ansatz' 

= e — (44^ 2 c) 1 (3 - 8) 

*i = e-^i / X ^**" W I (3 9) 
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with the two-spinors 3^.^ and ^J^- A general solution of the Dirac equation can be 
written as a linear combination of these wave functions (this is because one can obtain 
every combination of spherical harmonics in the four spinor components). 

In the regions where the i-variable is time-like, we choose the Dirac matrices again in 
the form (2.1), whereby the function S is now given by 

S(r) 



(3.10) 



According to ( |1.4[), the formula for the Dirac operator is obtained by inserting the Coulomb 
potential into (|2.5[), 



G 



id e 
Sdi~S 
d 



i) 



s' 







+ 1 7 — + «7^— 



dtp 



(3.11) 



The identity ( |3.3|) allows to rewrite the angular derivatives of the Dirac operator in terms 
of the operator K. If we substitute the ansatz' Q3.8| ),(|jO|) into the Dirac equation and 
apply the relations fl3.5|) , (|3 . 6[) , and (|3.7|) , we obtain the two-component Dirac equations 



S dr® jkuJ 



-1 

1 



(OJ 



1 
-1 



2j + l 
2r 



1 

1 



m 



jku 



(3.12) 



In the regions where the t-direction is space-like, we obtain the generalization of ( |2.8[) for 
the Dirac operator; namely 



G 



i d 
S di 
+ 1^8$ + i^8^ 



\Sdt r 



iS d r + S - + % - S' 
r 2 



(3.13) 



We again choose the ansatz' ( |3.8|) ,(3.9). This gives the two-component Dirac equations 



dr jku} 



-1 

1 



(CO 



&j>)—±% 



1 \ 2j + 1 



1 



2r 



+ i 



m 



*iL • ( 3 - 14 ) 



4 Non-Extreme Reissner-Nordstrom Background 



In this section, we consider the case q ^ p, so that the metric coefficient S(r), ( p. 10 ), has 
two zeros 



ro = P — y P q and n = p + y p 2 — q 2 ■ 

These zeros are transversal, S'(rj) ^ 0; in addition, the potential </>(r) is regular at r = rj. 
Since our matching conditions (|2.31 ) and ( 2.34 ) for the Schwarzschild metric only depend 



on the local behavior of the external field around the horizon, they are also valid for 
the Reissner-Nordstrom horizons (for the inner horizon, we must reverse the r-direction). 
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We will show in this section that these matching conditions do not admit normalizable, 
time-periodic solutions of the Dirac equation. More precisely, we will show that for every 
(non-trivial) solution of the Dirac equation ( |1.4[) , the normalization integral outside and 
away from the horizons, 

(*[*)^, := / S- 1 d 3 x , (4.1) 

Jm 3 \B 2ri 

is infinite for some t. Notice that for a normalized wave function, the integral (^^) gives 
the probability that the particle lies outside the ball of radius 2r\ , which must be smaller 
than one. Thus, if ( |4.1| ) is inifinite, the wave function cannot be normalized. 

Suppose that we have a periodic solution (|l.3|) of the Dirac equation with period T. 
Expanding the periodic function e *&(t, r, i?, ip) in a Fourier series gives the representation 
of ^ (as the Bloch wave) 

V(t,r,$,(p) = e- int *n(r, #,</>) e~ 2mn ^ . (4.2) 

n£Z 



Decomposing the functions ^f n in the basis ([T^), (pT9|) , and substituting into (4^2) gives 

*(t,r,#,<p) = ]T ^ikuwfrr^M > ( 4 - 3 ) 

n,j,k,s 

where the index s = ±, and where uj is related to n by 

. . 2im 
coin) = il + — — . 



Using the orthonormality of the two-spinors X ,i> the normalization integral takes the 
form 



x zr i n,n' j,fc,s 



t<Tis , % e-l ,73 

jku>(n') ° a x 



The integrand has an oscillating time dependence of the form exp(i(u(n) — u(n'))t). In 
order to eliminate the oscillations, we take the average over one period (0,T), giving 



n,j,k,s 

For a normalizable wave function ^, this expression is finite. Since the scalar product 
(•|-)oo is (semi-)positive definite, we conclude that all the summands must be finite; thus 

(*i*w(n) I *ifc«(n)) < 00 ( 4 - 4 ) 

for all s = ±, j, fe, n. 

This inequality allows us to turn our attention to the individual wave functions VP^,. 
As a first step we show that the wave functions $ in the ansatz' Q3.8| ) and ( |3. 9[ ) are not 
zero on the horizon. 

Lemma 4.1 T/ie function | < ^j^. w (^)| 2 ^os finite boundary values on the horizon. If it is 
zero on a horizon r = tq or r = r\, then & f kuJ vanishes identically. 
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Proof: For ease in notation, we omit the indices j, k, and oj. For a given 5, < 5 < vq, 
the i-direction is time-like in the regions (8, ro) and (7*1,00). In these regions, the Dirac 
equations ( |3.12 ) give 

si$±|2( r ) = < S V $±> + <$ ± ,^$ ± > 

dr dr dr 



1*11-1*2 1) " 4mRe(($f)*^ 
and thus 

_ c |$±|2 < 5 A|$±|2 < c |$±|2 

with c = 2m + (2j + Dividing by \^^\ 2 and integrating yields, for 5 < r < r' < ro, 

or n < r < r', the inequality 

S- 1 < logj^l 2 < c / S- 1 . (4.5) 



In the region ro < r < ri, the Dirac equations ( 3.14 ) give similarly 

5i$±| 2 (r) = <5^-$ ± , $±> + <$± 5^-d> ± > = , 
ar ar ar 

since the square bracket in ( |3.14[ ) is an anti-Hermitian matrix. Thus l^l 2 is constant in 
this region, and so, (|4.5| ) also (trivially) holds for ro < r < r' < r\. 

Notice that 5 _1 is integrable on the event horizons. Therefore, the inequality Q4.5| ) 
implies that the left and right sided boundary values of l^l 2 on the horizon are finite, and 
are non-zero unless "J? 1 * 1 vanishes identically in the corresponding region {8, ro), (ro,ri), or 

(n,oo). 

Next we consider the matching condition ( |2.34j ). If we substitute the ansatz' ( |3.§| ) and 
( |3 . 9|) , we get for the conditions 

$ ± (r i + e) - $ ± (r i - e) = o(l + |$ ± (r j + e)|) at e -> 0, j = 0, 1. 

Since we have already shown that ^^(r)! 2 has two-sided limits as r = r,-, this last equality 
shows that the left and right sided boundary values of l^^l 2 must coincide, 

lim \<& ± (r j + s)\ 2 = lim |$ ± (r 1 -e)| 2 ,7 = 0,1. 

We conclude that the wave function can only be zero on one of the horizons if it vanishes 
in the whole interval (<5, oo). Taking the limit 5^0 gives the result. ■ 

We point out that this lemma does not imply that the wave function <E> is continuous on 
the horizon. In general, $(r) will oscillate faster and faster as r approaches a horizon. 
Nevertheless, its absolute value |<3?| tends to a finite value in this limit. 

The next step is to use current conservation for analyzing the decay of ^j kLU ( n ) a * 
infinity. 

Theorem 4.2 (radial flux argument) Either ^ S j^ u vanishes identically, or the normal- 
ization condition is violated. 
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Proof: To simplify the notation, we again omit the indices s, j, k, and uj. Assume that 
\& is not identically zero. For n < r < R and T > 0, let V = (0, T) x (_B 2 ii \ B 2r ) be 
an annulus outside the horizon. As a consequence of the current conservation, the flux 
integral over the boundary of V is zero, thus 







v 

T 



alt r z Sir) / (^ 7 r ^)(t,r) 
J5 2 



dtR z S{R) / (^ 7 r ^)(t,i?) 



2fl 

'li- 



ds s 2 S-U 



s2 



(*7**)(t,r) 



t=0 



where / 52 denotes the integral over the angular variables. Since the integrand is static, 
the last integral vanishes, and we obtain that the radial flux is independent of the radius, 



r 2 S{r) / (^ 7 r *)(r) = R 2 S(R) / (tf 7 r tf)(ii) 



g 2 



s 2 



(4.6) 



We want to show that the radial flux is not zero. For this, we first substitute the 
ansatz' (ETH) and (|3,9p into the right side of (fO) and get 



r z S(r) / (^ 7 r $)(r) 
Js 2 



S 2 



i 
-i 



$(r) . 



(4.7) 



According to Lemma 4.1, |$| has finite, non-zero boundary values on the horizon ri. 
Expressed in <£, the matching condition ( p. 31 ) gives 



lim ( 



$ = . 



Using this equation, we take the limit r — ► n in (4.7), 



ri<r— >n 



lim r 2 5(r) / (^7 r ^)(r) 



s 2 



lim 



lim 

n<r—*ri J g2 



1 i 
-i 1 

1 
-1 



lim 

n <r^r 



/ |$| 2 ^ 
i Js 2 



<3? - |$r 



$ - i$i 



where we used Lemma [4.1| in the last inequality 

Now we consider the radial flux for large R. Since the flux is non-zero and independent 
of R, we have 

< lim \R 2 S(R) f (yyy)(R)\ . 



R^oo 



S 2 



Using the positivity of the form $7*^ and the fact that the Reissner-Nordstrom metric is 
asymptotically Minkowskian, we get (using the Cauchy-Schwarz inequality) the estimate 



< lim \R Z S{R) / (^7 r ^)(i?)| 



R->oc 



S 2 



R—>oo 



< lim \R 2 S(R) / (* 7 **)(i2)| = lim \R 2 S-\R) / (^7*^)(i?)| 



s 2 



R^oo 



S 2 
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We have shown that the integrand of our normalization integral 

roc r 

(*!*)«»= / dRR 2 S~\R) (W*){R) 

J2ri JS 2 

converges to a positive number. Thus the normalization integral must be infinite. ■ 

This theorem shows that the wave functions in the decomposition ( [4.2D , (|4.3| ) must 
all be identically zero. Thus there are no normalizable solutions of the Dirac equation; 
this proves Theorem |l.l| . 

Remark 4.3 We point out that the radial flux argument is based only on our matching 
conditions for the wave functions and on the Dirac current conservation. Therefore, it can 
immediately be applied to more general static, spherically symmetric background fields. 
This generalization may for example be relevant if the coupling of the gravitational and 
electric field to matter or other force fields is taken into account. Although the exact 
formulas of the Reissner-Nordstrom solution will then no longer be valid, the qualitative 
behavior of the fields on the horizons may still be the same. To give an example of the 
possible generalizations, we state the following theorem, which can be proved with very 
similar methods: Let gij be a static, radially symmetric background metric, 

ds 2 = g ij dx i dx j = ~^dt 2 - —— dr 2 - r 2 (d$ 2 + sin 2 d^ 2 ) , 



whereby the metric coefficient A{r) has N 
Assume the following conditions hold: 

(1) The zeros of A are all transversal, 



zeros at r = n, . . . , r^ , < ro < • • • < r^. 



A'( rj ) + forj = l,...,N. 

(2) The determinant of the metric is regular except at the origin, 

T- 2 {r)A~ l {r) £ C°°(0,oo) . 

Furthermore assume there is a spherically symmetric electric field cp(r) which is regular ex- 
cept at the origin, <f> G C°°(0, oo). Then there are no normalizable, time-periodic solutions 
of the Dirac equation with these background fields. 

5 Extreme Reissner-Nordstrom Background 

We now consider the case q = p of an extreme Reissner-Nordstrom background field, i.e. 

s= . 

r 

The metric coefficient 5 now has only one zero at r = p\ the £-direction is time-like both 
inside and outside the horizon. This situation can be thought of as the limiting case 
that the two horizons ro and r\ considered in the previous section come arbitrarily close. 
Unfortunately, the arguments for the non-existence proof do not carry over in this limit, 
so that we must rely on a different method. 
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Since the i-direction is always time- like, the i-component of the current ^G 1 ^ is pos- 
itive and has the usual interpretation as probability density. Therefore the normalization 
integral 



(* | *)* 



IR; 5 



causes no conceptual difficulties. 

Suppose that we had a normalizable, periodic solution (1.3) of the Dirac equation with 
period T. Again using the representation as the Bloch wave gll) and averaging over one 
period gives 

® n,j,k,s 

Substituting the ansatz' (|3.8|) and (|3.9|) yields 



' jkui(n) I jku>(n)> 



1 

T 



T 



(* I dt 



drS- 2 (r) £ |*$ M „)I 

n,j,k,s 



Using the positivity of the summands, we obtain the conditions 



dr S~ 2 (r) m 



jkiv(n) I 



< OO 



(5.1) 



for all s, j, k, and n. 

We will now study the individual functions for r > p. To simplify the notation, we 
again omit the indices j, k, and oj. Our first task is to consider under which conditions on 
the parameters u, j, and m the normalization integral (5.1) can be finite near r = p. We 
first discuss the situation qualitatively: Since S~ 2 (r) = r 2 /(r — p) 2 has a non-integrable 
singularity on the horizon, the normalization integral will only be finite if <3? s becomes 
small near r = p. For generic parameter values, the dominant term in the Dirac equation 
( |3,12| ) near r = p is the first summand, i.e. 



— *± 

dr 




-1 







4> d 



Since in this limiting case, the eigenvalues of the matrix on the right are purely imaginary, 
the Dirac equation describes fast oscillations of the wave function. The eigenvalues of the 
second and third summands in ( |3.12| ) are real; they describe an exponential increase or 
decay of If the oscillating term is dominant, we expect that $ will not go to zero in 
the limit r — > p. In the following lemma, these ideas are made mathematically precise in 
a slightly more general setting. 



Lemma 5.1 Let 3>(x), x > 0, be a nontrivial solution of the ODE 



&(x) 




-1 







+ b(x 








-1 



+ c(x 




$(x) (5.2) 



with smooth, real functions a,b,c € C°°(0, oo) and o^O. If near the origin, the quotients 
b/a and c/a are monotone and 

b(x) 2 + c(x) 2 < a(x) 2 , (5.3) 
then \$>\ 2 (x) is bounded from above and from below near x = 0, 

< liminf ^(x)) 2 < limsup [^(x)! 2 < oo . 



0<a;^0 
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Proof: Let (0, e) be an interval where the functions b/a and c/a are monotone and where 
( [5.3D holds. Assume that $ is a nontrivial solution of (|5.2|). According to the uniqueness 
theorem for the solutions of ODEs, <&(x) is non-zero for all < x < oo. Now consider the 
functional 



F(x) = <$(x), A(x) 4>(x)> 



with 



1 + 6/a —c/a 
—c/a 1 — 6/a 



According to ( |5.3| ), the matrix A is close to the identity; i.e., there is a constant c < 1 
with 

|1 — A(x)\ < c for all x with < x < e. 
Thus the functional F is uniformly bounded in |$| 2 on (0, e), 



(5.4) 



|$(:r)j 2 < F(s) < C |$(x)| 2 

for some C > 0. Using the special form of ^4 and of the differential equation ( |5.2D , the 
derivative of F takes the simple form 



F'(x) = <$', A<$>> + A$'> + v4'$> = A'$> . 
The sup-norm of the matrix A' is bounded by 

\A < 

Putting together (|5.4[) , (|5.5|), and (|5.6|) , we get the bounds 











+ 











(5.5) 



(5.6) 



a 



+ 



F(x) < F'(x) < C 



b\' 

a 



+ 



F(x) 



Now we divide by F(x) and integrate. Since b/a and c/a are monotone, we can just 
integrate inside the absolute values, 



C 





b 




c 




( 




+ 




) 




a 




a 





< log-F^ < C 





b 




c 




( 




+ 




) 




a 




a 





(5.7) 



Since the extreme left and right sides of this inequality converge in the limit x — > 0, we 
conclude that \ogF{x) is bounded from above and below near the origin. After exponen- 
tiating and substituting ( |5.4D , the result follows. ■ 



Applied to (|3.12|) , this lemma says that |$ ± (r)| 2 is bounded away from zero near r = p 
unless 

uj-e(j)(p) = . (5.8) 

Thus we can turn our attention to this special case. 

If we substitute the condition fl5.S{) into (3.12), the Dirac equation simplifies to 



1 



d_ 

dr 



$±( r ) 



-1 

1 



e ± 



1 \ 2j + l 
-1 j 2r 



1 

1 



m 



(5- 
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We want to study how the solutions of this equation behave for small r — p > 0. For this, 
we rewrite the equation in the new variable 



u(r) 



-r — p ln(r — p) 



which gives 



du 



$ ± (n) 



2j + 1 

2r 




(5.10) 



The region near r = p corresponds to large values of u. The matrix in the bracket in ( |5.10| ) 
depends smoothly on u and converges in the limit u — > oo to a finite limit , in view of the 



definition of u given above. According to the stable manifold theorem [14, Thm. 4.1], the 
solutions of ( 5.10) ) which are not bounded away from zero for large u tend exponentially 
to zero. After transforming back to the variable r, this justifies the power ansatz 



d?f(r) 



*To 



(r-p) s + o((r-p) s ) 



(r-p) s + o((r-pY) (5.11) 



with constants $^ , and a parameter s > 0. Substituting into (5.9) yields the system 
of linear equations 



(*=F(7 + 1/2))*m) 
(s ± (j + 1/2)) <D 



-p (m + e) $f Q 



± 

20 



-p (m 



e) $ 



± 

10 



(5.12) 
(5.13) 



which can be solved for and ^^q. In this way, we have found a consistent ansatz for 
the spinors near r = p. However, the corresponding solutions of the Dirac equation are all 
not normalizable, as the following theorem shows. 



Theorem 5.2 Every nontrivial solution $ ± (r), r > p, of the Dirac equation §5. 



with 



the boundary conditions (5.11) violates the normalization condition (5.1). 

Proof: Let be a nontrivial solution of the Dirac equation. Since the Dirac equation 
has real coefficients, we can assume that $ ± are real. In the new variable u = r , the 
Dirac equation (|5.9D takes the form 
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pu\ -^-<I> ± (m) 

du 








If e > m, Lemma 5.1 yields that ^^(tt)! 2 is bounded from above and below near u = 0. 
Thus ^^(r)! 2 does not decay at infinity, and the normalization integral (5.1) will diverge. 
We conclude that we must only consider the case m > e. 

In the case m = e, the syste m (p. 12 ), ( p.!3| ) yields that either $ 



or $20 ^ s zero - 



Furthermore, the Dirac equation (|5.9|) shows that either &f or c&if vanishes identically. 
Since ^^(r) has no zeros for finite r (otherwise, the uniqueness of the solution yields 
that vanishes identically), we can assume that the vector $^(r) will lie in the fourth 
quadrant, 

^(r) G {(x,y) | x > 0, y < 0} (5.14) 



for all r. 
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4>i 



/ ^ 




Figure 2: Flow of for large r, schematic 



Next we want to show that ( |5.14| ) also holds in the case m > e. In this case, from 
( j3p and (|5laD , we can assume that is positive, whereas is neg ative. Thus (15.14) 
holds for small r — p > 0. In order to show that the fourth quadrant is an invariant region 
for first notice that $ (r) cannot become zero for a finite value of r. Thus, if $ ± (r) 
leaves the quadrant for some r, we have either 



$f(r) = 



< and *jf (r) < 



or 



*f(r) > , $f=(r) = and (*a )'( r ) ^ 

But the Dirac equation gives in the first case that (^i 1 )' > and in the second case that 
(®f)'(r) < 0, which is a contradiction. 

We conclude that $ ± (r) lies for all r in the fourth quadrant. Figure [2] shows the flow 
of equation (5.9) for large r. From this one sees immediately that the origin is repelling, 
so that l^l 2 will be bounded away from zero for large r. ■ 

It follows that our periodic solution ^ must vanish identically outside the horizon. This 



proves Theorem 1.2. 



We point out that in contrast to the situation in Section U, we do not make any 
statement on the behavior of the wave function for r < p. Indeed, it appears that the 
extreme Reissner-Nordstrom background does admit periodic solutions for r < p; these 
can be constructed by taking the boundary conditions ( 5.111) on the horizon and solving 
the Dirac equation backwards in r. 
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A Justification of Time-Periodicity Inside the Horizon 



Throughout this paper, we have considered a Dirac wave function ( |1.3j ) which is time- 
periodic both inside and outside the event horizon. Since an outside observer has no 
knowledge about the physical situation in the interior of the event horizon, the assumption 
of time periodicity inside the horizon might not seem physically resonable. In this short 
appendix, we clarify why time periodicity inside the horizon is natural to assume. Namely, 
we show that every solution r,i?, <p) of the Dirac equation which is time-periodic 
outside the event horizon and (locally uniformly) bounded in t, gives rise to a solution ^> 
of the Dirac equation, which coincides with \& outside the horizon and is also time-periodic 
inside. Using this argument, the results of this paper could be immediately generalized to 
Dirac wave functions which are only time-periodic outside the event horizon. 

Let ^(t, r, if) be a solution of the Dirac equation which is time-periodic outside the 
event horizon, 

V(t + T,r,0,<p) = e~ inT ^(t,r,-d,(p) far r > n , (A.1) 

and locally uniformly bounded in t, 

\*(t,r,0,<p)\ < F(r) with FeC°((0,r )U(ro,ri)) (A.2) 

(ro and r\ again denote the Cauchy and event horizons, respectively). We consider for 
N > 1 the functions 

1 N 

* N (t,r,0,<p) = ^— Yl *(t + nT,r,>&,<p) . 

n=—N 

Since our Dirac operator is static, the functions ^> n satisfy the Dirac equation. Time- 
periodicity (A.l) implies that \£ n and coincide outside the event horizon. Inside the 



event horizon, one can use the bound ( |A.2| ) to show that the Vtjv form a Cauchy sequence. 
Thus we can take the limit N — > oo; we set = limjv-^oo ^N- Again using ( |A.2j ), we 
conclude that the function ^ is time periodic, 

V(t + T,r,0,<p) - *(t,r,0,p) 
= J im ^rrr W + (N + l)T,r,0,<p)-V(t-NT,r,0,<p)) = , 

N^co 2,1V + 1 

and satisfies the Dirac equation, 

(G-m)^> = lim (G - m) = . 
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